The perturbation theory of higher dimensional spacetimes 

d la Teukolsky 



Mahdi Godazgar 
Department of Applied Mathematics and Theoretical Physics 
Centre for Mathematical Sciences 
Wilberforce Road, Cambridge CBS OWA, UK 
2 ■ mmg31@cam.ac.uk 
(N 

_ October 27, 2011 

O 

o 

^ ■ Abstract 

(N . 

We consider the possibility of deriving a decoupled equation in terms of Weyl tensor com- 
^ ■ ponents for gravitational perturbations of the Schwarzschild-Tangherlini solution. We find a 

particular gauge invariant component of the Weyl tensor does decouple and argue that this cor- 
^ ■ responds to the vector modes of Ishibashi and Kodama. Also, we construct a Hertz potential 

map for solutions of the electromagnetic and gravitational perturbation equations of a higher 
dimensional Kundt background using the decoupled equation of Durkee and Reall. Motivated 
by recent work of Guica and Strominger, we use this to construct the asymptotic behaviour of 
^ ■ metric perturbations of the near-horizon geometry of the 5d cohomogeneity- 1 Myers-Perry black 

Q^ ' hole. 
I> ■ 

^ ! 1 Introduction 

Soon after the discovery, by Kerr, of a solution to the vacuum Einstein equation representing an 
isolated rotating black hole [Ij, the status of its classical stability spawned an area of active research 
in general relativity. Significant progress in this direction was made by Teukolsky, who realised that 
in the case of algebraically special solutions [21 E], of which the Kerr solution is an example, one is able 
to derive a decoupled equation, satisfied by Weyl scalar ^'oi from the original perturbation equation 
' [9l[10]. The validity of a separability ansatz, which is related to the existence of hidden symmetries 
[m [12] , allowed Press and Teukolsky to use the decoupled equation to provide strong evidence for 
the linear (mode) stability of the Kerr solution under non-algebraically special perturbations [13| . 

The existence of a decoupled equation was later shown to be related to the existence of gauge 
invariant quantities |14j . Given that the Weyl scalar ^'o, which solves the decoupled equation, is 
gauge invariant and has the same number of degrees of freedom as the perturbed metric, viz. 2, 
one may suspect that ^'o encodes all information regarding the perturbation, i.e. a solution of the 
perturbation equation can be constructed given the existence of a decoupled equation. This was 
shown to be true by Kegeles and Cohen [15[ [16] and Chrzanowski |17j . They outlined a constructive 
procedure (the Hertz map) for finding solutions of perturbation equations for an algebraically special 
background given some technical assumptions. A very short and elegant proof of these statements 
was provided by Wald some time later [18]. 

The study of higher dimensional gravity, particularly its black hole solutions, has become an 
active area of research in recent years. One of the most natural questions that one can ask is 
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whether known black holes solutions, such as the Myers-Perry [T9], or the black ring [20J solutions 
are classically stable. The stability of the Schwarzschild-Tangherlini solution, which can be thought 
of as belonging to the Myers- Perry family has been demonstrated |21j . 

However, it has been known for some time, by analogy with Gregory-Laflamme (G-L) type 
instabilities [22], that one would expect instabilities to occur in certain other regimes of the Myers- 
Perry and black ring families. For example, one would expect thin black rings to suffer from the same 
type of G-L instability that occurs for a black string [20]. Or, for ultraspinning (d > 6) Myers- Perry 
black holes to suffer from the same kind of G-L type instability that one finds for p-branes [23] . 
There has been a great deal of recent progress in using numerical methods to tackle and confirm 
these conjectures, at least, as far as the Myers- Perry solution is concerned [24j . 

Nevertheless, experience from 4d GR suggests that there may be a more simple framework (d la 
Teukolsky) in which the stability of higher dimensional black hole solutions can be addressed. The 
motivation of such a consideration would not only be that such a framework would facilitate a much 
simpler study of perturbations of regimes currently under investigation using numerical methods but 
that it may allow a study of regimes that are currently inaccessible to numerical investigations. 

The higher dimensional generalisation of the Teukolsky decoupled equation was first considered 
in Ref. [29]. In higher dimensions, one can construct the analogue of ^'o by choosing a null frame 
with null vectors i and n such that £ ■ n = 1 and complete this frame with d — 2 orthonormal 
spacelike vectors m(j) (i = 2, . . . ,d — 1) that are orthogonal to i and n. Now, we define the higher 
dimensional generalisation of ^'o to be iljj = ('^'m(^i^^£'^m(^j^'^Cabc^- Furthermore, for algebraically 
special solutions the perturbed value of Ojj is gauge invariant (under infinitesimal diffeomorphisms 
and frame transformations) ^29j. Thus, as the higher dimensional generalisation of \I'o, ^ij is the 
most natural candidate to consider decoupling in higher dimensions. 

By studying the conditions required for iljj to decouple, it was found that in higher dimensions, it 
is not sufficient for the background solution to be algebraically special. In addition, there must exist 
a geodesic null congruence with vanishing expansion, shear and twist. Equivalently, the solution 
must be Kundt. 

With the above motivation in mind, the strong restrictions found in Ref. [29], which a background 
must satisfy for ^lij to decouple, can be thought of as being disappointing. However, this result 
can still be used to study perturbations of a particular class of higher dimensional black holes. 
This is because the near-horizon geometry (NHG) of extreme black hole solutions are known to be 
Kundt. This fact was used by Durkee and Reall to find instabilities of the NHG of cohomogeneity- 
1 Myers- Perry solutions |30j . A conjecture relating instabilities of the NHG and the full solution 
for perturbations preserving certain symmetries was then used to predict an instability of extremal 
and near-extremal cohomogeneity-1 Myers- Perry solutions in seven or more dimensions [30] . These 
predictions were later confirmed in Ref. [31j . 

Even so, one would like a decoupled equation that allows the study of a greater class of solutions 
within the family of the Myers-Perry solution. While, it is true that Qij does not decouple it may be 
possible to construct some other gauge invariant quantity from the Weyl tensor that does. The first 
aim of this paper is to consider this possibility for the Schwarzschild-Tangherlini background as a 
simple class within the Myers- Perry family. The existence of gauge invariant quantities constructed 
from the Weyl tensor that decouple on the Schwarzschild-Tangherlini background would give us hope 
of finding analogous quantities for more general solutions within the Myers-Perry family. Equiva- 
lently, the absence of such quantities for the Schwarzschild-Tangherlini solution would indicate their 
absence for more general solutions. 

^Qij is the higher dimensional generahsation of "to in the sense that they are both the boost weight +2 components 
of the Weyl tensor. 
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In section O we begin by considering whether other gauge invariant quantities can be added to 
Ojj, so that the new quantity decouples. The strategy we take is to find the obstruction to the 
decouphng of fijj. Then, we consider the decouphng of this new quantity and find the obstruction 
to its decouphng. Once this iterative process terminates we are left with a set of gauge invariant 
quantities that form an obstruction to each others' decoupling. For the case of Qij, there are three 
basic gauge invariant quantities including fijj that must be considered. However, we find that a 
linear combination of these quantities can never decouple. Hence, there is no hope of constructing 
a gauge invariant quantity that decouples using Ojj. 

We then consider other gauge invariant quantities constructed from components of the Weyl 
tensor and find that a particular set of components <I>^ = £"n^m([j)'^m(j])'^Caf,cd do decouple. Metric 
perturbations of the Schwarzschild-Tangherlini solution have been studied by Ishibashi and Kodama 
|21j using a gauge invariant analysis that is analogous to the 4d gauge invariant approach developed 
by Moncrief [32]. This approach uses the spherical symmetry of the solution to construct gauge 
invariant quantities from the metric perturbations. These quantities can be classified into tensor, 
vector and scalar modes depending on how they behave on the sphere. We argue that is related 
to the vector modes of Ishibashi and Kodama. 

This is completely analogous to what is known to happen in 4d. On the Schwarzschild back- 
ground, the perturbed value of the imaginary part of the Weyl scalar ^2, which is gauge invariant 
and has $^ as its higher dimensional generalisation, satisfies a decoupled equation. Furthermore, 
this decoupled equation is equivalent to the Regge- Wheeler equation [34] describing vector mode (or 
axial) perturbations of the Schwarzschild background [33]. 

Our results decrease the likelihood of finding a gauge invariant quantity constructed from Weyl 
tensor components that decouples on the more general Myers-Perry background. This is because <I>^ 
is generally non-zero for the Myers-Perry solution, so its perturbation is not gauge invariant. It is 
not inconceivable that there exists a gauge invariant quantity constructed from that decouples. 
However, we have not been able to come up with a suitable candidate. 

The second aim of the paper is to construct the Hertz potential map for constructing solutions 
of the perturbation equations of higher dimensional Kundt backgrounds in the manner proposed by 
Wald [Hj. 

In section [3l we begin by presenting Wald's argument for constructing solutions of a general 
perturbation equation given a decoupled equation. It was shown in [29] that a decoupled equation 
for electromagnetic and gravitational perturbations exist for Kundt backgrounds. Thus, we apply 
this method to construct the Hertz map for electromagnetic (section 13. ip and gravitational (section 
13. 2p perturbations of Kundt backgrounds in higher dimensions. 

As stated above, the NHG of extremal black hole solutions is Kundt. In particular, the decoupled 
equation for NHG of cohomogeneity-1 Myers-Perry solutions was studied in [3^. As an application, 
we use the results of [3D] to determine the asymptotic behaviour of metric perturbations of the NHG 
of the 5d cohomogeneity-1 Myers-Perry black hole. This is, in part, motivated by a recent paper 
concerning the entropy counting of such black hole solutions [36]. We find that there exist modes 
that violate the boundary conditions required in Ref. [36j . Thus, at higher orders, it will not be 
possible to deform the NHG. 

2 Decoupling results for perturbations of the Schwarzschild- 
Tangherlini solution 

In this paper, we define the gauge invariance of a perturbed quantity following Stewart and Walker 
[14] . Given a spacetime scalar X, we can decompose this into its background and perturbed value. 
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i.e. to first order X = X*-*^-* + X^^\ We are interested in the how the perturbed value of X, X^^-* 
changes under infinitesimal coordinate and frame transformations. We say that X^^-* is gauge invari- 
ant if and only if it remains invariant under coordinate and frame transformations. This condition 
places constraints on the background value of X. For example, under a coordinate transformation 
parametrised by i, X^^) X(^) + ^ • dX^^\ Hence, for X(^) to remain invariant under infinitesimal 
coordinate transformations X^'^'^ must be constant. The invariance of the perturbed quantity X^^^ 
under infinitesimal frame transformations is verified by performing infinitesimal frame transforma- 
tions given in equations (|X2]) — (|X5]) of appendix E In Ref. [29], it was shown that the perturbed 
value of Ojj is gauge invariant if and only if the background solution is algebraically special. Similar 
arguments to those presented in Ref. [29] can be used to demonstrate the gauge invariance of other 
Weyl components given a background. For brevity, we will suppress superscript labels on quantities 
indicating whether they are background or perturbed values since this will be apparent from the 
context. 

In Ref. [29], Durkee and Reall consider the status of Teukolsky's decoupling result in higher 
dimensions. They do this by taking an algebraically special background, in which the perturbed 
value of rtjA is gauge invariant. Using the Bianchi identity equations of the higher dimensional GHP 
formalisiro [37], they derive a second order coupled differential equation, where all second order 
derivatives act on the gauge invariant quantity ^ij. For iljj to decouple, all other terms involving 
other components of the Weyl tensor must either be re-expressed in terms of terms involving Vtij or 
vanish. Once all the Newman-Penrose and Bianchi equations have been exhausted, they conclude 
that for VLij to decouple, the background must be Kundt. This means that the background admits 
a null geodesic congruence with vanishing expansion, shear and twist. 

Here, we shall take a different approach. Instead of forcing a particular gauge invariant quantity 
(in their case Vtij) to decouple and finding restrictions on the background, we fix the background 
(in this case the Schwarzschild-Tangherlini solution) and try to find gauge invariant quantities con- 
structed from Weyl tensor components that decouple on this background. 

Thus, we will use the higher dimensional GHP formalism developed in Ref. [37] to address the 
possibility that there exists a gauge invariant quantity constructed from Weyl tensor components 
that decouples on the Schwarzschild-Tangherlini background El- 

We begin by considering the derivation of the decoupling result in [29j for the Schwarzschild- 
Tangherlini background. Up to equation (5.31) in p9] . the only assumption made is that the 
background admits a null geodesic congruence with vanishing shear and rotation. This is satisfied 
by the Schwarzschild-Tangherlini solution. Removing other second order terms from equation (5.31) 
gives 

(2t)'l) + 6,6, + p'b + nP\^' + ^,PP' + + IS - ^,5,,) = 0. (2.1) 

Thus, the trace-free part of is the obstruction for the decoupling of iljj 0. We know that this 
vanishes on the background (see appendix [B]). Hence, it is invariant under coordinate transforma- 
tions. By applying infinitesimal frame transformations (IA.2P — (IA.5P and using the fact that the 
background solution is type D, we can simply verify that — -^^Sij is a gauge invariant quantity. 
We proceed by finding the obstruction to the decoupling of the gauge invariant quantity 

^ \ d-2"i-3 ) ■ 

^See appendix |X] for a review of the higher dimensional GIfP formalism. 

^See appendix |B] for a review of the Schwarzschild-Tangherhni solution. 

*Note that in 4d this term vanishes identically leading to the decoupling of ^o- 
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It can be shown that this quantity satisfied 

+ {S^P^P"^^^+P'^^,)=0■ (2-2) 

Hence, the obstructions to the decouphng of [^f, — jt^*^*? ) p^p''^^ij and p^^i^-, which can 



^ ij d-2 "V ^ 

simply be shown to be gauge invariant. An equation of the form above can be derived for these 
quantities using equation (j2.ip . However, this procedure will lead to yet another new quantity to 
consider: p^p'^ _ ^^-^-^ 

Thus, we find that this iterative process generates boost weight +2 terms of the form rijj, 
— and p'^^'ij multiplied by various factors of pp' . Therefore, we take 



X,, = n,, + fp' (^^>^^ - ^6,j) + g p'n',^ (2.3) 

as our ansatz for the decoupled gauge invariant quantity, where / and g are GHP scalars of boost 
and spin weight 0. 

For Xij to decouple, there must exist / and g such that 

{2])'l> + 6k6k)Xij = OXij, (2.4) 

where O is some first order differential operator. 

The only boost and spin weight GHP scalars that do not vanish on the background are <^ = pp' 
and Hence, 

/ = /(?,$), 5 = 5(?,$). (2.5) 

In particular, equations (|NP3p . (jNP3])' and (|B5P evaluated on the background imply that on the 
background 

6i/ = 0, 6,5 = 0. (2.6) 
Using equation (j2.ip . its primed version for p'^^'ij and eq. ()2.2p . we find 

(2i)'t) + 6fc6fc) Xi, = {p'\> + + 

+ { (-21^ + + f ) b + {-2^i +np)i>'+--] [fp' (4- - ^'^.•); 

+ { (-2i? + ^P' + f ) b + (-2lf + ^p) I.' + . . .} {gp'n[^) , (2.7) 

where the ellipses indicate other scalar terms that have been omitted for brevity. Necessarily, for 
Xij to decouple, the coefficients of the derivative operators in the 3 terms on the right hand side 
must be equal, i.e. 

P ~ ^ / d-2P "^p — g ^ d-2 P ^ p ' 

d-2P — f ^ d-2P — ^ s ^ d-2P- y'^-^) 

These imply that on the background 

b/ = 0, 1)5 = 0, ]3'/ = 2(sl2P'+f)/, 5 = 4 (h^p' + f ) 5. (2.9) 



^The derivation of this equation is given in appendix [Cl 
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Equation (lB2p evaluated on the background implies that 

J)'$ = -^p$. (2.10) 
This, along with the equations above, can be used to show that 

= {jj^PP' + '-^^ - ^) /, (2.11) 

which is non-vanishing on the background. However, considering commutator equation (jCip on /, 
on the background, gives that this should vanish. Thus, we encounter a contradiction. We would 
arrive also at a contradiction had we considered a commutator on g. 

Therefore, there do not exist scalars / and g such that the gauge invariant quantity Xij, given 
by eq. (|2.3p . decouples. 

In summary, we took as the initial ansatz for a decoupled equation. This does not work as 
there is another gauge invariant quantity that obstructs the decoupling of Ojj. Then, we considered 
the decoupling of this new quantity and found other quantities that obstruct its decoupling. Re- 
peating this iterative process we found that all obstructions are of a form given by one of three basic 
gauge invariant quantities constructed from ^fj and O^^. Thus, we take a gauge invariant quan- 
tity constructed from a linear combination of these three basic quantities as our new ansatz and find 
that such a quantity cannot decouple. Hence, a decoupled equation for the Schwarzschild-Tangherlini 
solution cannot be found by completing Qij with other gauge invariant quantities. 

It is clear from the above calculation that gauge invariant quantities formed from or the trace- 
free part of will not decouple. However, one can consider other gauge invariant quantities. For 
example, one can consider boost weight -|-1 components of the Weyl tensor: ^ijk- Although ^ijk is 
gauge invariant under infinitesimal coordinate transformations, since it vanishes on the background, 
its gauge invariance under infinitesimal tetrad transformations is not so clear. Under a null rotation 
about n parametrised by Zi, and ^ijk transform as |37) 

Zi<P + ^ijk^'i^ijk + Z[k^j]i + zi<^iijk + ---, (2.12) 

where the ellipses refer to terms that vanish on the background and need not be considered. Thus, 
in order to be gauge invariant, the boost weight -|-1 ansatz must be of the form 

^^ijk+d^Siy^k], (2.13) 

which vanishes in 4d. Looking at the Bianchi equations ()Bip - ()B7p . it is not clear how they can be 
manipulated to give a decoupling result for the quantity above. 
Alternatively, we can use the Penrose wave equation [3S] 

g'^VeVfRabcd + RabefKi + 2{RaecfRl{ " RaedfRl{) = (2.14) 

to confirm that decoupling is indeed not possible. Shortly after the decoupling result of Teukolsky, 
it was shown by Ryan, that the Teukolsky decoupled equation can be derived by projecting the 
Penrose wave equation on the appropriate null tetrad component and linearising [39] . Note, that 
since we are going to consider this equation to first order in perturbed quantities and, furthermore, 
that the perturbed solution continues to satisfy the vacuum Einstein equation, we can replace the 
Riemann tensors in the equation above with Weyl tensors. 

Contracting the Penrose wave equation ()2.14p with i"-m(^i-^^m(^j-^'^m(^f.-^'^ + ■^j^£""n''£'^5iym(^k]-^'^ gives 

o {^^Jk + = ^ {^'^}k - ihkjV) + • • • ' (2-15) 
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where we have used the Bianchi identities to simphfy the expression on the right hand side. The 
elhpses indicate terms that are of zero order in derivatives. The Bianchi identities can not be used 
to transform the first order in derivatives obstruction on the right hand side into a terms involving 
^ijk + d^^i[j^ k]- Thus, we conclude that this boost weight +1 gauge invariant quantity cannot 
decouple. 

We will not consider gauge invariant quantities constructed from ^ijki since it is not clear whether 
such a quantity would be simpler to work with than the original metric perturbation. 

The only components left to consider are It is simple to verify that is gauge invariant. 
We proceed by using the Bianchi equations to derive an equation, in which second order derivatives 
act only on <I>^. As before, we throw away terms that are clearly of quadratic order or above in the 
perturbation expansion. 

Antisymmetrising over ij in (1B2P gives 

- \><^tj + = s^P^ft + ^2M^Jh (2-16) 
while antisymmetrising over ij in (jBSP gives 

Contracting over indices jl in (IB3P and antisymmetrising over the remaining indices gives 

- ^k^m = - (2.18) 

Now, consider 8t. (f2l7D -2b' (l2T6D 

{21p'\) + bkbk)^tj + \>'^k^[ij]k - 2t)'6[i%] + 6[Ak\^j]k - ih', + = • • • , (2.19) 

where the ellipses refer to terms on the right hand side. All such terms have only single derivatives. 
In addition, commutator equations ()C2p and ()C3P can be used to rewrite the final two terms on the 
right hand side as an expression involving only single derivatives. Furthermore, equations ()C.4p and 
(|2T8]1 imply that 

- 2l3'6[i^'j] + 6[i6|fc|$j]fc = -[^',bii\^j] + ..., 

which can be converted into an expression with single derivatives using (IU2])'. Thus, we have an 
equation in which all double derivatives act on 

i2\)'\> + bkbk)<^>tj = ...■ 

On the right hand side, terms of the form 6fc^'[jj]fc can be removed by using equation (j2.18p . Other 
terms can be simplified by rearranging terms such that one of the factors in a term vanishes on 
the background. Then, the other factor can be assumed to take its background value. Bianchi and 
Newman-Penrose equations on the background can be used to further simply such terms. 

Having simplified terms on the right hand side of the equation above, we find that the expression 
simplifies greatly and decouples 

(2b'b + bkbk + + p't,) + ^^pp' - Hig) $A ^ 0, (2.20) 

Metric perturbations of the Schwarzschild-Tangherlini solution have been studied by Ishibashi 
and Kodama [21j using a gauge invariant analysis that is analogous to the 4d approach developed 
by Moncrief [32]. <I>^ is a 2- form on the {d — 2)— sphere. However, their analysis does not find 
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any 2-form type modes on the sphere. They classify the perturbations into three types of tensorial, 
vector and scalar modes depending on their behaviour on the (d — 2)— sphere, which is parametrised, 
here, by i type indices. The tensor modes decompose into a trace-free transverse symmetric tensor, 
a transverse vector and two other scalars, while the vector modes decompose into a transverse vector 
and a scalar. 

Performing a similar decomposition of gives 

^fj=^kj+h^j], (2-21) 

where ^^ij transverse, i.e. 

bi^^ij = (2.22) 

and Vi is some vector mode. Since transverse 2-form modes on the sphere are not found in the 
Ishibashi and Kodama analysis, we can only conclude that ^xij parametrise trivial metric pertur- 
bations. An intuitive reason for why this should be the case is that to construct (symmetric) metric 
perturbations from ^Tij^ ^^^^ need to contract it with a dj derivative, which gives zero. This 
leaves Vi, which correspond to the vector modes of Ishibashi and Kodama. 

As mentioned in the introduction, this is analogous to what is found in 4d. On the Schwarzschild 
background, the perturbed value of the imaginary part of the Weyl scalar ^'2, which is gauge invariant 
and has as its higher dimensional generalisation, satisfies a decoupled equation [33]. Furthermore, 
this decoupled equation is equivalent to the Regge- Wheeler equation |34j describing vector mode (or 
axial) perturbations of the Schwarzschild background [33]. The real part of the Weyl scalar 
which has $ = Cabcd as its higher dimensional generalisation, is not gauge invariant. Prom 

this, one can construct a gauge invariant quantity using metric perturbations and show that this also 
decouples and the decoupled equation is equivalent to the Zerilli equation ^35j describing scalar mode 
(or polar) perturbations of the Schwarzschild background |40) . Analogously, $ is not gauge invariant 
on a Schwarzschild- Tangherlini background. However, we do not attempt to construct a gauge 
invariant quantity from this since we are only considering gauge invariant quantities constructed 
from the Weyl tensor. Had we done so, we would presumably find that this satisfies the scalar mode 
master equation of Ishibashi and Kodama [21J. 

In conclusion, we find that the only gauge invariant quantities constructed from Weyl tensor 
components that decouple on the Schwarzschild background are ^fj. Purthermore, these components 
are shown to be related to the vector modes of Ishibashi and Kodama [21j. 



3 Constructing solutions of perturbation equations 

The existence of a map in 4d for constructing exact linear perturbations given the existence of a 
decoupled equation was demonstrated, given some technical assumptions, by Kegeles and Cohen 
[15\ [TB] and Chrzanowski [17]. However, Wald put this map on a firmer basis with his elegant proof 
of it [18]. In this section, we review Wald's result [18] regarding the construction of solutions of 
perturbation equations, given the existence of decoupled equations, emphasising the generality of 
his argument. 

We consider as our background, a solution of the vacuum Einstein equation 

Rab = ^-gab- (3.1) 
We would like to find a solution of some linear (perturbation) equation 

£if) = (3.2) 
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on this background, where / is the perturbation of some field. For example, for an electromagnetic 
perturbation, / is the 1-form potential A and the equation above reduces to the Maxwell equation 



V"(V„A- = 0. (3.3) 

In the case of gravitational perturbations, / is the associated metric perturbation, usually denoted 
by h, while £ is the linearised Einstein operator 



[0'^'Hh)U = ^cV^ahf - ^{V^Kb + VaVth) - ^gab (v,Vrf/l='^ - V^/l) + ^A{hgab - dhab) , (3.4) 

where all indices above have been raised using the background metric g, V is the Levi-Civita con- 
nection of g and h is the trace of hab- 

Now, assume that we can construct another equation, which is satisfied by a quantity T{f), 

O(r(/)) = 0, (3.5) 

where O and T are some linear operators. Such an equation, where T(/) is more simple than /, 
is usually referred to as a decoupled equation. In 4d, a decoupled equation for electromagnetic 
and gravitational perturbations can only be derived for vacuum solutions admitting a null geodesic, 
shear-free congruence [9ill0j. By the Goldberg-Sachs theorem [41j, this is equivalent to the condition 
that the solution be algebraically special (i.e. Petrov type II or more special). 

In higher dimensions, the existence of a decoupled equation for electromagnetic and gravitational 
perturbations places more restrictive conditions on the background solution [29j. For d > 4, the back- 
ground solution must admit a null geodesic congruence with vanishing optics, i.e. the background 
solution must be Kundt. Kundt solutions are necessarily of CMPP [32] type II |43] . 

The existence of a decoupled equation implies the existence of a linear operator S such that 

S£ = OT. (3.6) 

Simply put, operator S describes how the decoupled equation is derived from the perturbation 
equation. 

The key idea in Wald's derivation is the use of operator adjoints. Given a linear differential 
operator O acting on a tensor field t G T and taking it to a tensor field u £ U, the unique adjoint 
operator is defined via 

{u,Ot) = {0^u,t) (3.7) 

up to a total divergence, where the inner product between any tensor fields f i and V2 of rank m is 
defined by 

{V1,V2) ^ I Vi''^-^^V2a,...a^d{Vol). (3.8) 

We proceed by taking the adjoint of the operator equivalence (j3.6p and using the fact that 
(OiOs)^ = 0\0\ to get 

£^S^ = T^Ol (3.9) 
The operator equivalence above implies that given a solution ■0 of 

= 0, (3.10) 

S^ip solves the equation 

£\f) = 0. (3.11) 
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For all cases of interest in this paper, the operator £ is self-adjoint, i.e. f"^ = £. Thus, 

Ot^ = =^ S,(S^) = 0. (3.12) 

The solution ip is called the Hertz potential. Furthermore, 

0{TS^iP) = {Or)S^^ = {SS)SH = S8{S^^) = 0, (3.13) 

where we have used (j3.6p in the third equality and (|3.12p in the final equality. Hence, T5^^ solves 
the decoupled equation 

Oip = 0. (3.14) 
3.1 Electromagnetic perturbations of Kundt solutions 

In [29j, it was shown that for electromagnetic perturbations, a higher dimensional analogue of the 
decoupled equation found in 4d by Teukolsky [10] can only be derived for Kundt solutions. The 
higher dimensional decoupled equation is analogous to the Teukolsky equation in the sense that in 
both cases 7~(^) is the highest boost weight component of the Maxwell field F = dA. Hence, 

TE{A)i = ^^=t rat (VaA - VbAJ. (3.15) 

The decoupled equation is 

OE{^)i = 0, (3.16) 

where 

OE{v)i = {2\>'lp + bjbj + p'\> - Arjbj + (I. - 2^ a) + 2(-2r[,6,] + <i>| + 21>tj)v,. (3.17) 
As noted before, 

£E{A)h = V'^(VaA - VhAa) = 0. (3.18) 

Operator 5^; is found by considering the derivation of the decoupled equation in |29) . To simplify 
the derivation of Se, we assume from the beginning that the background under consideration is 
Kundt, i.e. = p^'^^j = 0. 

Also, in the derivation in [29] . a vector potential A is not introduced. Hence, the Bianchi identity 
dF = is non-trivial. Here we let F = dA. Thus, of the Maxwell equations (4.4)-(4.7) in [29J, only 
(4.4), (4.6) and their associated primed equations are non-trivial: 

(4.4) ^ £E{A)a = 0, (3.19) 
(4.4)' ^ -n'' £EiA)a = 0, (3.20) 
(4.6) ^ (4.6)' ^ m,'' £E{A)a = 0. (3.21) 

The derivation begins by considering the combination ]p(4.6) +bj [(5jj(4.4) — (4.5)], which results 
in equation (4.9) of [29j | 

= (2l3'i) + bjbj)ip^ + 2[l3, lp']ip^ - [b, 6,-] (F,, + F6^j) + [b,,bj]ipj 

+\>{ - (2/>(. .] - p'Sij)^j + 2{F,j + F5,j)tj) (3.22) 
-bi{T'jipj) +bj{2Tl.ipj]), 



®We use the notation of [29] to denote components of the Maxwell field strength; that is </?i = Pmi^i^ Fab, F 
tn^Fab, Fij = m^ifrrnj)''Fab and i/^i = n''m(^i)''Fab. 
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where, as noted before, we assume k = pij = 0. This is equivalent to considering 



(3.23) 



The next manipulation that involves equations (4.4) and (4.6) is eliminating the combination 
\){Fij + F6ij). This is done by adding -(2tj + T'j){6ij{4A) - (4.5)) to (|3:22]l or equivalently ([323]). 
Hence, we have 

l> (m,-^ £E{A)a) - bi {t £E{A)a) + {2n + t[) {t eE{A)a) . (3.24) 

The rest of the derivation does not make use of the Maxwell equation, that is it involves either 
the Bianchi identity or the Newman-Penrose equations. Hence, operator Se is given by equation 
([3:2i]l 

Se{J\ = l3(m,Va) - (6i - 2Ti - T'i){tJa). (3.25) 



The Hertz potential V'H is given by solving 

0^j^{^H)i = 0. 
We derive by considering the inner product 

{iI,,,Oe{v)^). 



(3.26) 



(3.27) 



For the inner product to be well-defined ^ipi must be boost weight -1 since OE{f)i is boost weight 1. 
Using equation ()3.17p 



ii^i,OE{v)i 



+ 2 



2d-3 
d-1 

S o^T>A 



A 



(3.28) 



Using equations (|C.3p - ()C.7p for the adjoint of the GHP covariant derivatives, commutator (|Cip and 
the NP equation ()NP4p gives 

OU^), = (21p'\> + bjbj + p't) + $ - ^) V'i + 2(-2r[,6,.] + ^ _ 2$A)^.. (3.29) 

To find a solution of the electromagnetic perturbation equation (j3.18p . we also need to calculate 
the adjoint of S. We do this by considering the inner product 



{^P^,S{A)i) = ([m,V+^"(-6^U27^, + T')' 
where we have used (j3.25p . Using equations ()C.3P and ()C.7p gives 

'S\lp)a = [-ma b + L{bi + Ti)] Ipi. 



i^i,Aa 



(3.30) 



(3.31) 



Thus, using the results of section [3l we have that if iI^h,^ is a solution to (j3.26p . where is given in 
equation (I3.29p . then 

sUi'H)a = [-m,a b + 4(6i + r,)] i^Hi (3.32) 

is a solution of the electromagnetic perturbation equation (jS.lSp . 
Also, 

TESl{ll}H)i = \^\^H)i (3.33) 
is a solution of the decoupled equation OEi^)i = 0. 
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Consider a doubly Kundt solution. This is a solution with two null geodesic congruences with 
vanishing optics, i.e. we also have k^^^^ = pf^^ = 0. In this case, the boost weight -1 component of 
the Maxwell field also satisfies a decoupled equation — the prime of the decoupled equation (|3.17p 

0'E{y^')i = 0, (3.34) 

where (using equation (jS.lTp and commutator (jCip ) 

O'Ei^'h = + 6,-6, - 2{tj + rj)6, + 2r,Tj - ^ - ^ a) ip', 

+2 [-2Tl^b,] + <^fj + + arf.r,]) (3.35) 

In this case, if tp'^j. is a solution to the prime of equation (j3.26p . then S'^\'ip'^)a is also a solution 
of the perturbation equation. 

3.2 Gravitational perturbations of Kundt solutions 

For gravitational perturbations we assume that the perturbed solution also satisfies the vacuum 
Einstein equation. The metric perturbation h satisfies 

ScWab = 0, (3.36) 

where Sq is given in equation (13.40 . 

In [29] , it was shown that as with electromagnetic perturbations an analogue of the Teukolsky 
decoupled equation for gravitational perturbations ^lOj only exists for Kundt solutions. The decou- 
pled equation is solved by the boost weight +2 components of the perturbed Weyl tensor which 
generalises the complex Weyl scalar ^0 of the 4d NP formalism to higher dimensions. 

TG{h),j = = -\ { (6(.6,) - ^6,8,) - 2 {r[fi,) - -f^r'.b,^ - (p^, - ^p') b 

+ - - (K(.)} {tm^'ha,) - ^ {\>bk - rl\> - {\>t'^] {tmk'Kb). (3.37) 
The decoupled equation is 

Oainh = 0, (3.38) 

where 

Oainhj = {^b'b + SfcSfc + P'l> - Qnbk + 4«> - ^a) Q^j 

+4 (Tfc6(,| - r(,|6fc + <^l^f^ + 4$^|fc)l7fc|^.) + 2<^ikji ftki. (3.39) 

To derive operator Sq, we assume that the gravitational perturbation also generates a first order 
energy-momentum tensor. Then, 5^ is the operator acting on the first order energy-momentum 
tensor in the inhomogeneous decoupled equation 

Ocmi, =87r Sg{t^^).^. (3.40) 

Therefore, we need to rederive the decoupling result of Ref. [2^ assuming a non-zero first order 
energy-momentum tensor. Thus, we must use the more general {viz. including matter) NP equations, 
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Bianchi identities and commutator relations found in Ref. [37J. In order, to simplify the calculation, 
we may assume from the onset that the background is Kundt. 

Going through the derivation of the decoupled equation for gravitational perturbations given in 
Ref. \29\, except using the more general equations that include matter terms gives 

Sg (T), . = ^ 5ij {2\>'\> + bkbk + p'\> - 6Tkbk + 4<^>) (n'Tab) - ['^l - \>p[^J)) {(''('Tab) 

+2 (i)6(, - (2r(, + T(',)t) - Opr'a)) {tm^'Tab) - t)l)(m,'^m/r,,) + ^ 6,, ma^'Tab). (3.41) 
The adjoints of Oq and So can be derived in a fashion analogous to the electromagnetic case 

ol{u),j = (2i3']3 + bkbk + p'l> + 2Tkbk + 4«> + ^^a) n,,- 

+4 (rfc8(i, - r(,|6fe + ^f.j^ - 4$(^|;- )nfcj^) + 2^ikji ^ki, (3.42) 

4 (n)a6 = -(aib ($1 + p[^j)lp) ^ij + 2^(a"^|i|f>) {\^^ + in + tI)\>) U,, - m,„m,-, ])^ U,, . (3.43) 
Thus, given a Hertz potential fin that is a solution of 

Olinnh = 0, (3.44) 
where operator Oq is given in eq. (I3.42p . then 

Sl{nH)ab (3.45) 

where Sq is given in eq. (I3.43p . is a solution of the gravitational perturbation equation p.36p . 
Furthermore, 

TgSI {^h\, = l^'^H^j (3.46) 

is a solution of the decoupled equation ()3.38p . 

In the case of a background that is doubly Kundt, the prime of all the above equations hold also. 

For the case of doubly Kundt solutions, the operators O, C, and O'^ encode the same physical 
information. Therefore, there should be some connection between them. Indeed, in 4d, for type D 
solutions, it is known that the solution of the equations corresponding to the above operators are 
related by various factors of the background value of the Weyl scalar ^2 (see Refs. |441 145} [18] and 
references therein). It is very simple to prove such relations using the Bianchi identities. However, 
in higher dimensions the situation is more complicated. We have not been able to derive similar 
expressions relating the solutions to the various equations defined by the four operators above. 

3.3 Asymptotic behaviour of metric perturbations of near horizon geometry of 
5d cohomogeneity-1 extreme Myers-Perry solutions 

In this section, we consider as an application of the Hertz potential map developed above, the 
asymptotic behaviour of the metric perturbation of the near horizon geometry of 5d cohomogeneity-1 
extreme Myers-Perry black hole. This consideration is motivated by a recent proposal that quantum 
gravity on the near horizon of a class of 5d solutions of which the above solution is an example, with 
appropriate asymptotic fall-off conditions on the metric perturbation is equivalent to a CFT, which 
can be used to calculate the Bekenstein-Hawking entropy of the original solution [36]. Thus, giving 
a statistical counting of the black hole's degrees of freedom. 
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Turning off tlie graviplioton charge in the solution discussed in [36] gives the 5d cohomogeneity-1 
extreme Myers-Perry black hole solution with near-horizon geometry 

ds'^ = -R^dT'^ + ~R2- + + cos 9d(t) + RdTf + dn^^^ \ , (3.47) 

where r+ is the event horizon radius and d^^^g) round metric on S"^. This (doubly Kundt) 

solution is studied in Ref. [30], where the decoupling result of p9] is used to predict an instability 
of the corresponding extreme Myers-Perry solution. Choose the following null frame 

-RdT + —], n = RdT + — 



1712 = ^[d'4) + COS e dcj) + RdT), rua = — e^, (3.48) 
v 2 2 

where a = 3, 4 and form an orthonormal basis on S"^. 

In Ref. [30], it was shown that for the geometry with metric (j3.47p with basis chosen as above 



Ki = K'i = 0, pij = p'ij = 0, Tj + = 0, 

4 

^ij = n'ij = 0, ^ij-fc = ^ijfc = 0, <^^jkl = RijkU ^Tj = 2" ("^3 A mi)ij, (3.49) 

T _i 



where Rijki is the Riemann tensor of the three dimensional space % with metric 



„2 



ds^ = ^ |2(dV' + cos ed(t)f + c?J^(2)} • (3.50) 

In order to determine the asymptotic behaviour of the metric perturbation of the near horizon 
geometry with metric (I3.47p . we must first solve the Hertz potential equation 

Oj.(n)y = 0, (3.51) 

where operator o'q is given in (I3.42p . Then, 

Kb = ^ iah{RijUij) + 2i(^am\j\b)lp6iUij - miarrij^ J)^ Uij (3.52) 

is a solution of the gravitational perturbation equation in the ingoing radiation gauge, where we 

have used eq. (j3.43p and eqs. (|3.49p above. Space % is not an Einstein solution. Thus, the la^b 
component of hab is non-zero. 

Assume the following separability ansatz for Iljj 

Iiij = x{T,R)Y,j{d,<i),i:). (3.53) 

Substituting this ansatz into the equation for Iljj, i.e. eq. (j3.44p gives 

{D^ -q^ -X)X = ^, 0(2)Ky. = AFi,-, (3.54) 



where D is some charge covariant derivative for some AdS2 scalar with charge q defined in Ref. [30] 
and O^^-* is some operator given in Ref. [30], which we do not need to know about in detail here. 
Hence, x solves the equation for a massive, charged scalar in an AdS2 background with homogeneous 



'^The steps involved in this computation are almost identical to that given in appendix A of Ref. [30] 
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electric field. Such an equation has been studied by a number of authors [Ml STJ SU [30]. At large 



(3.55) 



Hence, from the form of hab given in (j3.52p we can conclude that for large R 

foil) 

<^iw) ^(tf) (^iw) 



hab ~ ^2 



2^ 2'l 



\ 



Vl + 4A. 



(3.56) 



where the columns and rows indicate the T, R, ip and 9 and </> (collectively labelled a) components, 
respectively. The same result would be found if we considered n as the WAND of choice, that is if 
hab were in outgoing gauge. 

Comparing this with the fall-off conditions in Ref. [36], we find that the TR, Tip, RR, ipa and 
a/3 components are the most restrictive. These components have been underlined in the matrix 
above. Thus, in order to satisfy the fall-off conditions, r/ must be real. We must also choose the 
lower sign (corresponding to normalisable modes) and require that > 1 or A > 0. Recall that A is 
the eigenvalue of operator 0(2)_ xhe spectrum of operator O^^) was studied in appendix B of Ref. 
[30j . It is clear from the study of gravitational scalar modes that there exists modes for which 



X = 2 + k{k, + 1) + \m\{K + \)- 



(3.57) 



where m is an integer and k is a positive integer. For k = 1, m > 15, A < 0. Hence, there exist 
modes that violate the fall-off conditions. It was shown in Ref. [30] that for all axisymmetric modes 
(m = 0) A > 0. Hence, all axisymmetric modes satisfy the boundary conditions. 

There is a similar proposal for the entropy counting of the 4d extremal Kerr solution [l9]. In 
Ref. [49J, as in the 5d case, the metric perturbation of the NHG is assumed to satisfy a certain 
asymptotic form. The asymptotic behaviour of the NHG of the extremal Kerr solution has been 
studied in [57] and [35]. The results we find for the 5d case are the same as they found in the 4d 
case. 
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A Higher dimensional GHP formalism 



In this appendix, we review the higher dimensional GHP formahsm of [37j. Given a background 
solution, we choose a null frame (£, n,m(j)) such that in this frame, the metric takes the form 

gab = 2^(„nb) + m(i)^m(i)^. (A.l) 

In the GHP formalism, one breaks complete covariance by singling out two null directions {I and 
n) at each point, but preserves covariance in the remaining directions. This is in contrast to the NP 
formalism where none of the covariance is preserved. 

At any point, the Lorentz group divides into 

• boosts (A a real function): 

^— >-A^, n — )• A^'-'-n, m(j) — t- m(j) , (A. 2) 

• spins [Xij £ SO{d-2)): 

i ^ £, n ^ n, m(^i) — Xijm(^j-), (A. 3) 

• null rotations about i {zi d — 2 real functions): 

n — )• n + Zim(^i) — \z'^£, m(^i) — m(j) — Zii, (A.4) 

• null rotations about n {zi d — 2 real functions): 

^ — )• ^ + ZiTn^^i-^ — ^z'^n, n ^ n, m(j) — )• — ZiU, (A. 5) 

where A 7^ and Xij is some position-dependent orthogonal matrix. 

We would like to keep the subgroup that preserves the null directions, i.e. the subgroup given 
by boosts and spatial rotations, or spins. Thus, we would like to work with objects that transform 
covariantly under this subgroup. 

Define any scalar rji-^...i^ that transforms covariantly as 

r]n...is >^^rih...is (A.6) 

under boosts and 

under spins, a GHP scalar of boost weight b and spin s. Evidently, the product of two GHP scalars 
of boost weights 61 and 62 and spins si and S2, respectively, gives a GHP scalar of boost weight 
61 + 62 and spin si + S2- 

Defining the covariant derivatives of the basis vectors as 

i 

Lab = Vfe4, Nab = ^bna, Mab = Vbm^i)^, (A.8) 

one finds that not all the scalars formed from the projection of these objects into the basis are GHP 
scalars. Those that are GHP scalars are listed in tabled] [37j . 

Notice that we have used a prime operation, which interchanges the null basis vectors 

' : i^n. (A.9) 
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Table 1: GHP scalars constructed from covariant derivatives of the basis vectors. 



Spin coefficient 


GHP notation 


Boost weight h 


Spin s 


Interpretation 






1 


2 


expansion, shear and twist of t 


La 


n = Da A 


1 





p"y"DansinTi of / 






2 


1 


non-geodesity of I 


Lii 


n 





1 


transport of I along n 




p'ij 


-1 


2 


expansion, shear and twist of n 


Nii 


P' = P'^^ 


-1 





expansion of n 


Nil 




-2 


1 


non-geodesity of n 




< 





1 


transport of n along / 



The prime operation is especially useful when considering type D backgrounds, since in this case I 
and n are essentially equivalent. 

The non-GHP covariant scalars formed from the covariant derivative of the basis vectors can be 
used to construct GHP covariant derivatives . For a GHP scalar rji-^,„is of boost weight b and spin 
s, we define its GHP covariant derivatives to be[f| 



\>Ti^i2...is = ^ ■ f^^iii2... is ~ ^-^10^4112. ..is + Afv0^ii...i,— ifcir+i... is' (A. 10) 

r=l 

Til 42... is = ^ • dTi-^i^ i^ — hLiiTi-^i^ i^ + MviTii...i^_ifcir+i...»si (^-H) 

r=l 

6j7jij2...is = ^{i) • ^^iij2...is ~ bLuTjj^j^^^^j^ + M>i^ji...>-ifcjr+i-.-js- (^-l^) 



In GHP notation, the Newman-Penrose, Bianchi and the commutator equations are much more 
compact that in the NP formalism [50]. For convenience, we write these equations, here, for an 
Einstein spacetime |29j . 

A.l Newman-Penrose equations 

\)Pij-bjKi = -pikPkj - KiTj - TiKj - Qij, (NPl) 
\>Ti-l>'Ki = p^ji-Tj+T!j)-^i, (NP2) 

26y|Pi|fe] = 2TiPijk] + '^i^iP'ijk] - ^ijk, (NP3) 

\)'pij-bjTi = -TiTj - KiKj - PikPkj - ^ij - -^Sij. (NP4) 

Another four equations can be obtained by taking the prime ' of these four (i.e. by exchanging the 
vectors i and n). 



Symbols J) and 8, pronounced "thorn" and "eth", respectively are old Germanic letters that have been retained in 
the Icelandic alphabet. 
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A. 2 Bianchi equations 
Boost weight +2: 

-2i^[jl6ii + ^i6[jii + + ^[j\ii)pi\k] + 2J^i[jr^], (Bl) 

Boost weight +1: 

+{^ijk - ^iSjk)T'k - 2{^^i6j)k + ^(ii)fc)rfe - ^ikp'kj, (B2) 

+4^^tjP[kl] - 2<^[k\iPj\l] + 2^[k\jPi\l] + 2^ij[k\mPm\l] 

-2^imT{j] - 2^ik\ijr'i] - 2niik\p'jii^ + 2n^ikP-\i], (B3) 

-^[j\^i\kl] = 2$y^|Pj|i] - 2^iyPkl] + ^im[jk\Pm\l] " 2ilj[jPfc;], (B4) 

Boost weight 0: 

\>'^ijk - 2bij\<^i\k] = 2{^-^6ii - ^j\u)Pi\k] + {2^i[jSk]i - 2Sii^ff, - ^iijk)Ti 

+2{^,6[jii - ^^lJ\l)p^k] + 2f^»[j4]> (B5) 

-26[i«'ffe] = 2^'|,p,fc] + |/5;|fc]-2^'[,p;.,]-^'i[,,|/5;|,], (bg) 

-^[k\'^ij\lm] = -"^'ilkllPjlm] + %[kl\Pi\m] " 2'^[k\ijP\lm] 

-"^iMPjlm] + - 2*[fe|iiP|i„^]- (B7) 

Another five equations are obtained by applying the prime operator to equations (IBip - (IB5p above. 
A. 3 Commutators of derivatives 

Acting on a GHP scalar of boost weight b and spin s, commutators of GHP derivatives can be 
simplified by: 



s 

+ - + W - ^v^i + 2^tj) (CI) 



r=l 



[t, 6i]Tki...ks = -{l^i\>' + T-b + Pjibj)Tk^...ks + b {-'T'jPji + l^jP'ji + ^i) Tkr...ks 

s 

+ J2i'^k-P'li- Pkrin +^krPli- P'kri'^l-^ilkr)Tkl...l...ks, (C2) 



r=l 



[bi,^j]Tk,...ks = (2p[,^]b' + 24]i)) Tk,...k^ + 6 (2pz[.|p;|,] + 2$^) Tfc^...^^ 

+ (^^'^rWP/lj] + '^p'kr[i\Pl\j] + ^ijW + Tk^...l...ks- (C3) 

r=l 

The result for [])' , bi] can be obtained from dCS])'. 
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B Schwarzschild-Tangherlini solution 

The Schwarzschild-Tangherlini black hole is an example of a higher dimensional type D solution. 
The Schwarzschild-Tangherlini metric in Schwarzschild coordinates is 

ds^ = -fir)df + dry fir) + r^dnf^_,), f{r) = 1 - (B.l) 
where dd'^^_2^ is the round metric on a unit radius (d — 2)— sphere. The WANDs of the solution are 

m 

5r ' dv ~^ ^"^ dr ^ ^ 

Defining (d — 2) orthonormal spacelike vectors 

TUi = rei (i = 2, • • • ,d), (B.3) 

completes a null frame {£,n,mi) for the solution. 

Cartan's first equation of structure de^ + u^^ye" = can be used to find the optical scalars 
associated with WANDs H. and n 

Lii = —Noi = —^drf, 
Ki = Ti = 0, pij = -j^Sij, p= [d- 2)/r, 

4 = rl = 0, p',^ = ^6,,, p' = {d-2)f/{2r). (B.4) 

Thus, the solution is an example of a Robinson- Trautman solution, that is there exists a null geodesic 
congruence with vanishing shear and rotation, but non-vanishing expansion. 

The curvature tensors can be derived from Cartan's second equation of structure i?^,^ = doj^y + 
ijOf^f A ujpu. Or, alternatively, one can read off the curvature tensors from appendix A of [52] 

^ij = ^ijk = 0, ^'ij = = 0, 

^^^■'^'^ (d-2)(d-3)' ^^.--2^^^J^' % = 0' ^ = • 

In particular, 

^% - ^'^u- = 0. (B.6) 



C Derivation of equation (12.21 ) 

This appendix is dedicated to the derivation of equation (j2.2p . 

We proceed by deriving an equation in which second order derivatives act only on — ^r^'^y- 
To simplify the derivation we shall neglect from the beginning any terms that are clearly of quadratic 
order or above when quantities are decomposed into background plus perturbation parts. 

Contracting ()B3p . taking its symmetric part and removing its trace gives 

2]3 ($1 - + 8fc - - ^5,j) =p^kl) {^^kjl + S'^u) - P {'^il - 
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Symmetrising over the ij indices in (jBSp and removing its trace gives 



6^ ( ^1 - A'J^. 



P 
d~2 



(ti-l)<E> Z' ^ _ Tfc X 



(C.2) 



Now consider jj^ fjC.ip + bh <\C.2\i . the left hand side of which is equal to 



6fc ( 



d-2 "^J 



(C.3) 



The top line is precisely what we want since commutator equation ([C2])' can be used to convert the 
second term to an expression with only first order derivatives. The second line can be simplified 
using the trace of (IB5P and full contraction of ()B7p . which when added together give 



Applying 6j to this, symmetrising over ij and removing the trace gives 



(C.4) 



(C.5) 



Commutator equations ()C3P and ([C2])' can be used to rewrite the second line in terms of first order 
derivative terms. Thus, the equation above can be used to rewrite the second line of (IC.3P in terms 
of first order derivative terms. To summarise, we have 

{2\>'], + bubu){^%--^5,: 



(C.6) 

where the right hand side of the equation above involves only first order derivatives. 

Terms of the form 6fc'I'(jj)fc and bk^'^^-^j^ can be removed by using equation (|C.ip . while a term 
of the form \>'Q.ij can be removed by using the symmetrisation of ()B2p . 

The remaining terms can be simplified by rewriting them such that one has a factor that vanishes 
on the background. Since, we are neglecting terms of order quadratic or above, this means the other 
factor takes its background value. For example, a term of the form 



P{ki 



P{kl) - dZI 
2 



Ski ) V ( '^ikji + ^Sij ) + ^Ski\>' ( ^ik,i + 



"(d-2)(d-3) 



Pi^3)-dh^^J)^''^-i^vm-^,6. 



2<I>5 



2{d-l) 



(C.7) 



where in the second equality, we have used the fact that — -g^Ski vanishes on the background 
and since we are only considering terms that are of linear order in the perturbation expansion, we 
can take the background value of the factor multiplying this term. In the final line we use Bianchi 
equation (jB2])' evaluated on the background to simplify Using this trick to simplify all such 
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terms and, also, using (1NP4P to eliminate terms of the form 6jTj, we find that the equation simplifies 
significantly 

Of course, we would like to derive an equation satisfied by ~ Letting <I>^*j = 

p'^'\>^% = {p%^%) - {\^'pV\, = Kb {p'K) - {\^'p')]^^% - {\>P'W^% - {\^'\>P')^%, 



Equations (jNPip and (|NP2p evaluated on the background give 

\>P = -^,, Kp=-^-^. (C.9) 
Thus, equation (jC.SP is equivalent to equation (j2.2p 

(2t'i> + 6.6. + ^(p'b + + f b + - i^^i^H^'^) [p' - A^.); 

D Adjoints of GHP covariant derivatives 



In this appendix, we derive the adjoints of the GHP covariant derivatives. First, consider the adjoint 
of operator \>. Let r/j^ . j^ and Cii...is be GHP scalars of spin s and boost weights h and —{b + 1), 
respectively and consider the inner product 

-6i/loC«ii2...is ~^11t=1 Cii...ir-ikir+i...is i^ ' 9 Tji-^i^ i, 

s ^ 

-[{I - d +V ■ I) Ciyi2...is + bLloCi-^i2„,is + Y.r = l MirO Cil...ir-lkir + l...is],Vili2-i, 

k 

i...ir~ikir+i...is]i 'niii2---ii 
~ ~l" P]Ciii2---is ) Viii2---is) ) 

where the first equality uses the definition of operator given in eq. (jA.lOp . the second equality 
uses the property that 

Mjt, + Mi^ = (C.l) 

and the third inequality is obtained using integration by parts and ignoring divergence terms, since 
operator adjoints are defined up to such terms. The penultimate equality uses the geodesity of i to 
deduce that 

V-£ = Lio + p (C.2) 
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and the final equality uses the definition of operator \) given in eq. (jA.lOp . Hence, 

b^ = -(b + p)- (C.3) 
Taking the prime of this equation gives the adjoint of Ja' 

\>'^ = -{V + p')- (C.4) 

Now, consider the inner product of i)i^r]i^...is+i with ^i^.i^^j, a GHP scalar with boost weight —b 
and spin s + 1 



\...is+i Xyr=l -^«rn^n--V-lfcv+l---«s+l ' '^(ii) ■ ^^i2---«s+l 

■["^(n) • 9 + V • m(i^) + 6Liij - Mkk] ^n.-Wi 

~X/r- = l -^«rnCn--ir-lfc*r + l---«s + l! Vi2---is + l 

-[8ii +V-m(ij) - Mfcfc]^ii...j,+i,??j2...Wi^ 

= — Til ~ ''"iil'^n-.-is+n ^i2---i3+i) > 

where the first equality uses the definition of operator 6 given in eq. ()A.12p . the second equality 
uses the property that 

Mj> + Mi^ = (C.5) 

and the third inequality is obtained using integration by parts. The penultimate equality uses the 
definition of operator 8 and the final equality uses the fact that 

i 

V • m(j) = Mkk -Ti-T^. (C.6) 

Thus, 

= -6i + Ti + t[. (C.7) 
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